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ABSTRACT 



Aims. We aim to investigate the chemistry and gas phase abundance of HNCO and the variation of the HNCO/CS abundance ratio as 
a diagnostic of the physics and chemistry in regions of massive star formation. 

Methods. A numerical-chemical model has been developed which self-consistently follows the chemical evolution of a hot core. 
The model comprises of two distinct stages. The first stage follows the isothermal, modified free-fall collapse of a molecular dark 
cloud. This is immediately followed by an increase in temperature which represents the switch on of a central massive star and the 
subsequent evolution of the chemistry in a hot, dense gas cloud (the hot core). During the collapse phase, gas species are allowed to 
accrete on to grain surfaces where they can participate in further reactions. During the hot core phase surface species thermally desorb 
back in to the ambient gas and further chemical evolution takes place. For comparison, the chemical network was also used to model 
a simple dark cloud and photodissociation regions. 

Results. Our investigation reveals that HNCO is inefficiently formed when only gas-phase formation pathways are considered in 
the chemical network with reaction rates consistent with existing laboratory data. This is particularly true at low temperatures but 
also in regions with temperatures up to ~200 K. Using currently measured gas phase reaction rates, obtaining the observed HNCO 
abundances requires its formation on grain surfaces - similar to other 'hot core' species such as CH3OH. However our model shows 
that the gas phase HNCO in hot cores is not a simple direct product of the evaporation of grain mantles. We also show that the 
HNCO/CS abundance ratio varies as a function of time in hot cores and can match the range of values observed. This ratio is not 
unambiguously related to the ambient UV field as been suggested - our results are inconsistent with the hypothesis of Martin et al 
(2008). In addition, our results show that this ratio is extremely sensitive to the initial sulphur abundance. We find that the ratio grows 
monotonically with time with an absolute value which scales approximately linearly with the S abundance at early times. 

Key words. Astrochemistry - Stars: formation - ISM: hot molecular cores 



1. Introduction 

The molecule HNCO (isocyanic acid) was first detected in the 
interstellar medium over 30 years ago jSnvder & Buhl|[T972l) in 
the Sgr B2 molecular cloud complex where its distribution was 
found to be spatially extended with relatively strong emission. 
Since its discovery, HNCO has been subsequently observed in 
severa l other molecular clouds, including t he dark cloud TMC - 
1 (e.g. lBrownlll98llGoldsmith et al.lfTgSirjackson et al.l[T984l) 



In the survey of 18 molecular clouds conducted by Jackson et alj 
([T98l . HNCO emission was detected in seven sources with 
an average excitation temperatu re of 12 K, siiggestin g that this 
molecule was tracing cold gas. I Jackson et alj d 19841) proposed 
that HNCO was a dense gas tracer due to the coincidence of 
the emission with regions of high density (n > 10^ cm"^). 
A la rger survey of 81 molecular clouds was conducted by 
IZinchenko et al. (2000). They reported a 70% detection rate with 
typical fractional abundances relative to the total hydrogen num- 
ber density, hh, of ~ 10"^ and a wide range of rotational temper- 
atures from 10 K to 500 K. Three of the objects were mapped and 
showed the emission region to be compact and centrally peaked. 

Isocyanic acid has also been detected in hot cores 
(HCs), the sites of ear ly stage, high mass star-formation 



including G34.3-h0.15 (iMacDonald et all Il996h . W3(H20) 
(lHeImich& vanDishoecklll997h and Sgr B2: in the N,M and 



NW cores (i Nummefin etall l2000f). A survey to wards seven 
high-mass star forming regions by Bisschop et alj (|2007) iden- 
tified HNCO with fractional abundances relative to «h of be- 
tween 6.4x10"** and 5.4x10"^ and excitation temperatures from 
64 K to 278 K. The observations suggested that HNCO can 
also trace much warmer gas - consistent with the findings of 
Zinchenko et al. (2000) and as such could be a potential indi- 
cator of star formation activity. 

A recent multi transition study of 13 molecula r clouds to- 
wards the Galactic Centre was conducted by Martin et al.l 
dlOOl). The sources were selected to represent a range of phys 
ical conditions and included photodissociation regions (PDRs 
along with HCs and giant molecular clouds. Martin e t al.l (1200 
compared the CS abundance, a species used as a dense gas 
tracer and is regarde d as a tracer of star formation sites 
dBronfman et al J 19961) . to that of HNCO for each of the sources. 
They found that the ratio of the abundance of HNCO to that 
of CS (as traced by '■'CS) was smaller for regions in which 
the FUV radiation is enhanced (PDRs), perhaps indicating that 
HNCO is more sensitive than CS to FUV photodissociation. 
The HNCO/CS ratio was systematically smaller in PDR-like 
regions by approximately two orders of magnitude than in the 
giant clouds. HCs were found to have intermediate values. 
It was therefore proposed that the HNCO/CS abundance ra- 
tio may provide a useful tool in distinguishing between shock 
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and radiation (FUV ) activity in molecular clouds. However, the 
iMartm et all (l2008 l) survey only observed a limited sample to- 
wards the Galactic Centre and before the usefulness of this pos- 
sible tracer can be established more complete surveys and a the- 
oretical understanding of the processes involved are desirable. 

Since the detection of HNCO there has been much specu- 
lation as to how this molecule forms but there have been few 
chemical models developed which incorporate an HNCO reac- 
tion network. So me authors have considered forma tion only in 
the gas phase (e.g. [Turner et al.|[T999Hlglesiasll 1977b whilst oth- 
ers have modelled the molecule's formation on grain surfaces 
re.g. lGarrod et al.l2 008'). Here we reconsider th ese processes and 
motivated by the observations of iMartm et al.l (|2008). construct 
a model for a hot core including both gas phase and grain surface 
formation and destruction of HNCO. We set out to model both 
the HNCO abundance and the HNCO/CS ratio as a function of 
time for such a region. We also aim to explore the response of 
the gas-phase formation of HNCO to changing physical condi- 
tions. Besides modelling HNCO in hot cores, simple models of 
a dark cloud and PDRs are also presented for comparison. 

2. Model descriptions 

Our hot core model consists of two distinct stages. The first stage 
is the collapsing cloud (CC) phase which follows the isothermal 
collapse of a dark molecular cloud whilst the second stage fol- 
lows the evolution at constant density and a specified tempera- 
ture - this is the HC phase. Along with the CC and HC models 
simple dark cloud (DC) and PDR models were also produced. 
These models were similar to the HC models in that no collapse 
was included but different physical conditions were used. All the 
models described in this paper were single point models. 



2. 1. Physical model 

The CC model begins at a time t - with an homogeneous dark 
molecular cloud at a temperature, T and an initial total hydro- 
gen number density of 2x10"* cm"^. The total hydrogen number 
density is defined as «h = 2n(H2) + n(H) and is equal to the 
cloud's density. It was assumed that the cloud was shielded by a 
visual extinction. Ay, of 15 magnitudes and the standard Galactic 
UV radiation field and cosmic ray ionisation rates were adopted. 
These physical conditions were chosen to be representative o f 
a typical Galactic molecular cloud (e.g. 'Millar & Herbstl ll990h . 
The cloud was then allowed to collapse isothermally following 
the modified free-fall collapse formula (see lSpitzeiill978h : 



dt 



4\l/3 



lAnGm-nrio 



n(t) 



1/3 



- 1 



1/2 



cm"^s"',(l) 



where B i s the retardation factor and a value of 0.7 was adopted 
following iRawlings et al.l (Il992l) . «□ is the initial cloud density 
(total hydrogen number density) prior to the collapse being ini- 
tiated, G is the Universal Gravitational constant and niu is the 
atomic mass of hydrogen. 

The collapse was halted once rtn - 2 x 10^ cm"^ which oc- 
curs after approximately 5x10^ years. This density was chosen 
to be representative of a typical hot core. As soon as this density 
was reached the model terminated and this was considered to be 
the end of the CC stage. The final chemical states of these col- 
lapse models were used as the initial chemistry of the subsequent 
hot core phases. A summary of the physical parameters adopted 
for the seven CC models used in this work is given in Table [T] 



Table 1. Summary of the seven collapsing cloud models. T is 
the isothermal temperature of the collapse. 

Model T (K) HNCO formation chemistry 

CCl 10 All 

CC2 10 Gas phase only 

CC3 10 Surface only (reactions a,b) 

CC4 10 All, reduced sulphur 

CCS 20 All 

CC6 50 All 

CC7 100 All 



The HC model represents the second stage of the overall evo- 
lution. The HC was considered to be an homogeneous cloud 
at the final density reached at the end of the collapse phase, 
«H = 2 X 10^ cm"^, in which the Ay was high enough that the 
influence of external FUV photons is negligible. For the duration 
of the HC model these physical conditions remained fixed. 

A total of 10 HC models were produced. Models 1 through 
4 all had the same physical conditions including a temperature 
of 200 K, meant to represent the heating by a newly formed cen- 
tral star The four models differ in their chemical networks. On 
the other hand, models 5 to 10 all had the same chemistry but 
differing temperatures (the temperature of the initial the collaps- 
ing clouds for these models also differed). A summary of these 
models can be seen in Table |2] 

In order to follow chemistry on dust surfaces it was neces- 
sary to include dust grains in the models. By assuming a gas 
mass to dust mass ratio of 100, the dust number density (rela- 
tive to «h) was calculated to be 1.33x10"'^. A single population 
of spherical dust grains was then assumed with each grain hav- 
ing a radius of 10"^' cm. The total number of binding sites on 
each grain was approximately 10^ assuming a surface density of 
7.9xl0i'*cm-l 

The DC and PDR models were modelled in a similar fashion 
to the HC. The physical conditions were kept constant in each 
case and no collapse was included. The DC was modelled using 
a r = 10 K, «H = 2 x 10** cm"^ and an Ay = 15, whilst the 
PDR model adopted conditions representative of the Horsehead 
nebula (Habart et al. 2005); T = 50 K, hh = 2 x 10^ cm-\ Ay = 
1 and an enhanced FUV radiation field equal to 60 times the in- 
terstellar value as used in all the other models. In order to mimic 
the effects of self-shielding of H2 and CO, which is not explic- 
itly included in the code, the photodissociation rates for both 
these species were set to zero in the PDR models. All the DC 
and PDR models produced in this work are also summarised in 
Table|2] A range of PDR models were run to demonstrate the ef- 
fect of different aspects of the HNCO chemistry, different densi- 
ties and different initial chemical compositions, either atomic or 
evolved from a collapsing cloud. A very dense PDR with an ini- 
tially evolved chemistry such as models PDR4 and PDR5 were 
explored to represent clumps exposed to the UV radiation from 
a very nearby young massive star. 

2.2. Chemical model 

For the CC and DC models, along with PDRl, PDR2 and PDR3, 
all species were initially present in their elemental (atomic) 
form. The only exception was hydrogen which was assumed to 
be predominantly molecular with only a small fraction present 
as atoms. These initial abundances (relative to rtu) are given in 
Table [3 For the HC models, along with PDR4 and PDR5, an 
evolved initial chemistry was used. The initial abundances at 
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Table 2. Summary of the dark cloud (DC), hot core (HC) and PDR models used to model the time dependent HNCO gas-phase 
abundance. T is the gas-phase temperature and «h is the total hydrogen number density (both held constant). The initial chemistry 
was either atomic, where all species (except H) were in their elemental form, or evolved, where abundances from a collapsed cloud 
(CC) model were used. The collapsing cloud models are listed in a separate table. 



Model T (K) iih (cm Initial chemistry HNCO formation chemistry 



DCl 


10 


2.0x10^ 


Atomic 


Reaction c only 


DC2 


10 


2.0x10'' 


Atomic 


Reactions c,d,e and f 


DC3 


10 


2.0x10* 


Atomic 


Reactions d,e,f and g 


DC4 


10 


2.0x10* 


Atomic 


Reactions c,d,e,f and g 


HCl 


200 


2.0x10' 


CCl 


All 


HC2 


200 


2.0x10' 


CC2 


Gas phase only 


HC3 


200 


2.0x10' 


CC3 


Surface only (reactions a,b) 


HC4 


200 


2.0x10' 


CC4 


All 


HC5 


20 


2.0x10' 


CCS 


AH 


HC6 


20 


2.0x10' 


CCl 


AH 


HC7 


50 


2.0x10' 


CC6 


AH 


HC8 


50 


2.0x10' 


CCl 


AH 


HC9 


100 


2.0x10' 


CC7 


AH 


HCIO 


100 


2.0x10' 


CCl 


AH 


PDRl 


50 


2.0x10* 


Atomic 


AH 


PDR2 


50 


2.0x10* 


Atomic 


Gas phase only 


PDR3 


50 


2.0x10* 


Atomic 


Surface only (reactions a,b) 


PDR4 


50 


2.0x10' 


CCl 


All 


PDR5 


50 


2.4x10' 


CCl 


All 



the start of these models used the abundances taken from a CC 
model. All the HC models and PDR4 used abundances taken 
from a CC at the final density, 2 x lO' cm""*, whereas PDR5 
used abundances from a CC model at an earlier stage in the col- 
lapse. In this case the density reached in the collapse was equal 
2.4 x 10^ cm"^ 

For all the models, once the chemistry had been initialised, 
the subsequent evolution of all species was followed via a chem- 
ical reaction network. This network included both reactions in 
the gas phase and on dust grain surfaces, along with the coupling 
mechanisms of accretion and thermal desorption only. 

The gas-phase chemistry included all the reactions from 
the late st UMIST database for astrochemistry (UDfAB, 
Rate06 (IWoodall et al.1 l2007h . For the grain-surface chem- 
istry over 200 reactions were inclu ded from Allen & RobinsonI 
(ll977h . lHasegawa et al.1 (Il992h and lHasegawa & Herbst (Il993l) . 
Reactions on the surface were treated using the modified rate 
equation method as used by Ruffle & Herbst (2000). 

In order to more completely model the HNCO gas-phase 
abundance some additional chemistry was included. Five surface 
reactions involving HNCO were taken from th e expanded OSU 
gas-grain reaction network (iGarrod et al ]|2008) comprising of a 
single formation path (reaction b in Table |4]i plus four destruc- 
tion mechanisms. This surface chemistry was complemented 
by some gas phase reactions (also from iGarrod et al.l l2008l) 
which followed the destruction of four new species (HNCHO, 
HNCOCHO, HNCONH and HNCOOH) formed from surface 
processing of HNCO. The destruction of HNCO in the gas phase 
was al so included. This followed the treatment of iTurner et al.l 
(II99 9I) and included the reactions between HNCO and the fol- 
lowing species; H+, H+, He+, C+, H3O+, CO+, HCO+, OCN+ 
and HNCO^. The rate coefficients for these reactions we re taken 
from the ion-neutral reaction database compiled by AnicichI 
(120041) . The photodestruction of HNCO by FUV photons and 
cosmic ray i nduced photodissociation was also incorporated 
dGredel et al .' 1989; Rober ge et al..,1991i) . Reactions with ions. 



particularly Hj and He^ dominate the destruction of HNCO in 
both dark clouds and hot cores. 

Two gas-phase formation reactions for HNCO were added 
(reactions c and d in Table |4]l. The first reaction (c) is a neutral- 
neutral reaction whilst the second (d) is the dissociative recom- 
bination of HaNCO^ ■ For re action c the rate coefficient is poorly 
known. iTurner et al] d 19991) adopted a value of 10 cm^ s ' in 
their dark cloud model. However the energy barrier for this re- 
action may be much larger than suggested by this number and as 
such this reaction may only become significant in regions much 
warmer than dark clouds. The rate coefficient for this reaction 
has been measured but only at temperatures above 300 K (see 
the NIST chemical kinetics databaso- A fit to the available data 
provides an expression for the temperature dependent rate coef- 
ficient; 



k(T) 



1.5 x 10"" xexp 



-4465 \ 



3 -1 
cm s , 



(2) 



where T is the temperature. At 10 K this gives a con sider- 
ably smaller value than that used by [Turner et al.l (Il999l) . We 
used model DCl to examine how the rate of reaction c affects 
the HNCO abundance - the model was run several times, each 
time using a different value for the rate coefficient ranging from 
10"^^ cm-' s"' to 10"** cm"* s"' (increasing the value by two or- 
ders of magnitude each time). For all other models in which re- 
action c was included the rate was calculated using Equation|2] 
Three further reactions were added that form 'parent' species 
needed for the formation of HNCO. Reactions e and / are ion- 
neutral reactions that are the first two reactions in a three stage 
process. The actual formation of HNCO in this reaction chan- 
nel is via reaction d (Ilglesiaslfl977h . The rates of these reac- 
tions are calculated using the Langevin rate. The final reaction 
ig) is a charge exchange reaction between He^ and OCN which 
produces OCN^. This reaction represents a third, new, reaction 



UDfA; www . udfa . net 



^ The NIST chemical kinetics 

http : / /kinetics . nist . gov/kinetics/index . j sp 



database; 
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Table 3. Standard initial chemical abundances relative to the to- 
tal hydrogen number density; F(X) - n(X)/nu. 





Y(X) 




nx) 




K(X) 


Hz 


0.5 


H 


1.00 X 10-" 


He 


0.1 


C 


1.32 X 10-" 


N 


7.50 X 10-5 


O 


3.19 X 10-" 


S 


1.86 X 10-5 


Si+ 


2.46 X 10-^ 


Fe+ 


1.50 X 10-** 



Table 4. Key reactions included in the chemical network relevant 
to the formation of HNCO as included in the chemical network. 
Thejjrefix G denotes surface species. AR77:Allen & Robinson 

19771) . GWH08:Garrodetalj (l2008h . TTH99; iTurner et afl 

19991) . I77: llglesia&(.1977.) 



Label 


Reaction 


Reference 


a 


GOCN + GH 




GHNCO 


AR77 


b 


GNH + GCO 




GHNCO 


GWH08 


c 


OCN + H2 




HNCO + H 


TTH99 


d 


H,NCO+ + e- 




HNCO + H 


177 


e 


H2 + OCN+ 




HNCO+ + H 


177 


f 


H2 + HNCO+ 


—> 


H2NCO+ + H 


177 


g 


He+ + OCN 


— > 


He + OCN+ 





channel between these two species since the UDfA already con- 
tains two reactions between He^ and OCN. A rate coefficient 
of IQ-^ cm^ S-' was adopted for reaction g - representative of 
charge exchange reactions. 

All models were allowed to evolve chemically for 10^ years, 
except for the collapsing cloud models which were stopped once 
the final density was reached. 

3. Results 

3.1. Gas phase formation of HNCO in dark clouds 

For the first part of our investigation, using DCl, we explored 
the effect the rate coefficient for reaction c has on the abundance 
of HNCO in the gas phase. Figure [1] shows the HNCO abun- 
dance along with that of OCN for a range assumed rate coeffi- 
cients for reaction c. For a rate coefficient of IQ-'^^cm-' s-' and 
above it is possible to achieve HNCO ab undances within the ob- 
served range of 1-4 x 10" '° to 5.4 x IQ-^ jMacDonald et al.ll996t 
iBisschop etsU 'lOOT). Equation |2]imphes that for the ~ 4500 K 
barrier which laboratory data (NIST) infer for this reaction, tem- 
peratures above 277 K are required to reach this value of the rate 
coefficient for reaction c. Alternatively for this reaction to be ef- 
fective at a gas temperature of 10 K, the barrier would have to be 
less than ~ 165 K, much lower than measurements indicate. In 
other words, at the temperatures of dark molecular clouds these 
reactions do not occur 

Model DC2 illustrates how the other gas-phase formation 
route (via reaction d) influences the HNCO abundance. Reaction 
c was also included in this model and the rate was varied as 
before. Figure |2] shows the HNCO and OCN abundances along 
with those for OCN+. A lower limit for the fractional HNCO 
abundance (at steady state) is effectively set at 1.1x10-". This 
corresponds to the DC2 model with a rate cofficient of for re- 
action c (the solid red line in Figure |2]). In this model HNCO 
is underabundant compared to observations by over 2 orders of 
magnitude. However, as seen for DCl, when the rate coefficient 
for reaction c increases, it is possible to achieve higher abun- 
dances. There is however very little difference between DCl and 
DC2 for those models with non-zero rate coefficients for reac- 
tion c. Figure [3] shows the gas-phase abundance for the key re- 
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Fig. 1. Time dependent, gas phase relative abundances for 
HNCO ( top 

) and OCN (bottom) as a function of rate coefficient for reaction c in 
a non-collapsing dark cloud (model DCl). All other HNCO formation 
pathways are prohibited, both in the gas-phase and on grain surfaces. 
Reaction g has rate coefficient of 0.0. The shaded region defines the 
range of values for the observed HNCO abundance, 1.4 x 10-'" to 5.4 x 
10-^ 



actants/products involved in the gas-phase production of HNCO 
(model DC3). This shows that there is OCN in the gas-phase 
but due to the low temperature this is not converted into HNCO. 
Instead, at low temperatures HNCO formation through reaction 
d (via reactions e and /) clearly dominates. The fourth and final 
DC model, DC4, includes all the gas phase formation reactions 
listed in Table |4] (but in common with the other DC models, 
no grain surface processes) and the reaction rate for reaction c 
has been varied. The abundance profiles for this model are seen 
in Figure |4] For small values of the rate coefficient of reaction 
c, the inclusion of reaction g increases the HNCO abundance. 
However, to achieve HNCO abundances in the observable range 
for ages less than 10^ years, the Turner rate must be greater than 
10-'^ cm-' s-', which would imply temperatures of > 370 K. At 
ages greater than -10^ years reaction g appears to be able to 
produce sufficient HNCO to match the observations even for a 
much reduced rate for reaction c (indeed even a rate of zero). 
However interpreting this late time rise in HNCO abundance 
needs to be treated with caution as these models only include 
gas phase reactions and so do not include the effect on the gas 
phase abundance of freeze-out on to the grains which will be 
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Fig. 2. Time dependent, gas phase relative abundances for 
HNCO {top), OCN {middle) and OCN+ {bottom) as a function 
of rate coefficient for reaction c in a non-collapsing dark cloud 
(model DC2). All gas-phase HNCO formation pathways are al- 
lowed; grain-surface reactions are prohibited. Reaction g has rate 
coefficient of 0.0. The shaded region defines the range of values 
for the observed HNCO abundance 
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Fig. 3. Time dependent, gas phase relative abundances for vari- 
ous species in a non-collapsing dark cloud using only reactions 
d,e and / to produce HNCO (model DC3). 



3.2. HNCO in ptiotodissociation regions 

Figure |5] shows the results for our PDR models. The HNCO 
abundance in PDR regions are even lower than those found in 
DC. For the three PDR models with an initially atomic chem- 
istry (PDRl, PDR2 and PDR3) the largest late time abundance 
of HNCO (10 is produced in the models with gas phase for- 
mation (PDRl and PDR2), but these produce up to 9 orders of 
magnitude less HNCO than in a DC model. In model PDR3, 
which allows only surface reactions, HNCO is even less abun- 
dant, reaching a maximum value of 10"^". The final two PDR 
models are also presented in Figure |5] Both models follow a 
similar trend, that is the HNCO abundance is steadily decreas- 
ing with time, reaching a plateau at later times. However, the 
HNCO abundance in PDR4 is initially larger by about two or- 
ders of magnitude than in PDR5, although they ultimately reach 
similar abundances of about 2x10 '** . None of the PDR models 
ever approach the abundances of HNCO seen in DCs and HCs 
- even for the PDR models with start with previously evolved 
chemisti-ies (PDR4 and PDR5). 

The warmer conditions in these PDR models, largely due to 
the reduced optical depth and enhanced radiation field resulting 
in the heating of dust grains, lead to faster desorption meaning 
that once formed HNCO is never on the grains for long (at the 
50 K used in these models the surface lifetime is approximately 
100 years) so surface formation is limited. Furthermore, HNCO 
is efficiently destroyed by both FUV and cosmic ray induced 
photodestruction reactions. At no time during the evolution of 
any of our PDR models d oes the HNCO abun dance come close 
to the ~ 10 observed by lMartm etaP ('2008') towards PDR re- 
gions. As a consequence of this low HNCO abundance, the ratio 
of the abundance of HNCO to the CS abundance in these PDR 
models is al s o man y orders of magnitude lower than seen by 
iMartm et al.1 (l2008l) (F ig.lTll. The origin of the HNCO observed 
by Martin et al. ( 20081) apparently in PDRs is therefore unclear 



very significant at these late times. This peak in HNCO abun- 
dance might occur earlier in models of denser clouds, but the 
effect of freeze-out on the gas phase abundances would also oc- 
cur earlier 



3.3. HNCO in iiot cores 

Figure |6] shows the HNCO gas-phase abundance as a function 
of time for all the hot core models 1 to 10. Only times between 
1000 years to 1 Myr (since the HCs switch-on) are plotted as 
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Fig. 4. Time dependent, gas phase relative abundances for 
HNCO (top), OCN {middle) and OCN+ {bottom) as a function 
of rate coefficient for reaction c in a non-collapsing dark cloud 
(model DC4). All gas-phase HNCO formation pathways are al- 
lowed but grain-surface reactions are prohibited. Reaction g has 
rate coefficient of 1.0x10"^ cm^ s"'. The shaded region defines 
the range of values for the observed HNCO abundance 
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Fig. 5. Time dependent, relative gas phase abundance of HNCO 
for the five PDR models. For comparison the observed value is 
between 1.4 x 10"'° abd 5.4 x 10"'^. Note that curves for PDRl 
and PDR2 are coincident. 
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Fig. 6. Relative HNCO gas-phase abundance as a function of 
time for hot core models (HC) 1-4 {top) and 5-10 {bottom). 



HCs are transient objects and are relatively short-lived phases of 
high-mass star formation lasting less than 1 Myrs. 

HCl includes the complete reaction network described ear- 
lier The HNCO abundance for this model does reached observed 
levels. Beyond ~ 10"* yrs values greater than 10"' are seen and 



a peak of ~ 2x10"^ is reached at around 0.2-0.3 Myr For the 
model with surface formation only (HCS) the lower observed 
limit for HNCO is reached but only beyond 10^ yr. The HNCO 
abundance in HC4 (which differs from HCl only in the elemen- 
tal sulphur abundance) is almost identical in behaviour as HCl. 
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The major differences are that higher abundances are reached 
at shghtly earher times and the peak value of just over 2x10"*^ 
remains for a further 0. 1 Myr. 

The gas-phase formation only HC model (HC2) shows that 
the abundance for HNCO remains relatively flat and never ex- 
ceeds a few 10 and even then only at HC ages greater than 
0.3 Myr As discussed above, this indicates that even the rela- 
tively high temperatures in hot cores, the gas-phase reactions by 
themselves are incapable of producing HNCO at observed lev- 
els. 

The remaining HC models all had the same chemical net- 
works (the full network with all reactions 'on'), it was only 
the initial and final temperatures that differed. In fact for the 
sake of comparison three of the models, HC5, HC7 and HC9, 
each had a constant temperature during the initial collapse and 
throughout subsequent constant density evolution. The first two 
models, HC5 and HC6, had the same post-collapse temperature 
(20 K), however HC6 used evolved abundances from a cloud 
collapsed at 10 K wheras HC5 used abundances from a cloud at 
20 K. HC5 never has a fractional HNCO abundance above 10"'^, 
whilst HC6 never really exceeds 10 HC6 has consistently 
more HNCO throughout - due to a lower temperature collapse 
phase temperature. At lower temperatures the rate of thermal 
desoprtion is reduced such that the surface chemistry becomes 
more efficient. This means there is more NH, CO and OCN on 
grain surfaces - the key reactants for the surface formation of 
HNCO. For the models presented in this paper, the surface re- 
action between OCN and H is more important for the formation 
of HNCO than the reaction between CO and NH. 

The next two models had post-collapse temperatures of 50 K. 
Again, the only difference being the collapse temperature (HC7 
with r = 50 K, HC8 with T = 10 K). Both of these mod- 
els have higher HNCO abundances compared to HC5 and HC6. 
Model HC8 is almost exactly the same as HCl, but HC7 has a 
relatively constant fractional abundance of around 10"'^. Model 
HC9 had a collapse and post-collapse temperature of 100 K. In 
this model the abundance of HNCO is roughly constant. At early 
times abundances just greater than 10 " are seen, falling slowly 
to just above 10 beyond 0.1 Myr The final model, HCIO, also 
has a post-collapse temperature of 100 K but the collapse tem- 
perature is 10 K. This has an almost identical abundance profile 
as HC8. These results indicate that in order to achieve HNCO 
gas-phase abundances similar to those observed, post-collapse 
temperatures greater than 50 K are required. 

3.4. The abundance of CS 

Our models can also be used to investigate the abundance of 
CS as function of time in hot cores (and PDRs). A detailed dis- 
cussion of the important pathways for the evol ution of sulphur- 
bearin g species, including CS, is given in iMiUar & Herbi 
(Il990l) . As mentioned above, CS is often used as a tracer of 
dense gas associated with star formation an d lMartmetaf] (120081) 
suggested the use of the ratio of the abundance of HNCO to that 
of CS as measure of the importance of UV radiation in a region. 

Figures |2] and [8] show the evolution of the abundance of CS 
in our PDR and HC models respectively. All three PDR models 
which start with atomic compositions (PDRl, PDR2 and PDR3) 
show identical behaviour. From its initial value of 2 x 10"'** the 
abundance of CS monotomically decreases by a factor of > 10 
over the period covered by the model. For PDR4 and PDR5 
which start with evolved chemistries the CS abundance also 
monotomically decreases from a peak at early times, becoming 
approximately constant at later times. The highest abundance of 
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Fig. 7. Time dependent, relative gas phase abundance of CS for 
the five PDR models. 
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Fig. 8. Time dependent, relative gas phase abundance of CS for 
four hot core models (HCl-4). 



CS in these models occurs are very early times in PDR5, where 
the peak abundance of ~ 10"^. In the HC models the CS abun- 
dance is initially relatively constant, rising to a peak after about 
10^ yr before rapidly declining thereafter. Comparing models 
HCl and HC4, which has a sulphur abundance ten times smaller 
that the other models, shows that the CS abundance throughout 
most the evolution of the HC scales nearly directly with the ini- 
tial sulphur abundance. 

4. Discussion 

Hot core models HCl, HC2 and HC3 demonstrate that gas 
phase reactions alone can not produce a high enough abundance 
of HNCO to match observations. Reactions on the surfaces of 
grains are required. Of the models with a complete HNCO chem- 
istry only models HCl, 8 and 10 achieve observed HNCO abun- 
dances. (Model HC4 can also match the observed abundance, but 
only differs from HCl in its initial sulphur adundance.) These 
models all have a low collapse temperature combined with high 
post-collapse (hot core) temperatures. It appears that during the 
collapse phase cooler temperatures are essential. (A HC model 
identical to HC4, but with an inital collapse phase at 20 K, 
rather than 10 K, produces a peak HNCO abundance of only 
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~ 3 X 10 '^.) 1000 years into the HC phase, reaching a peak 
abundance at ~ 10^ years of only ~ 10 The cold collapse 
allows accretion on to the surface of grains and hence surface 
production of HNCO. However, during the hot core phase, high 
temperatures are needed to allow the desorption of the surface 
species back into the gas-phase. The models with constant tem- 
peratures across both stages 1 and 2 of the hot core model are 
seemingly incapable of producing HNCO at observed levels. 

Interestingly, although grain surface reactions are clearly the 
ultimate source of the HNCO, the HNCO abundance does not 
follow the expected time dependent behaviour of a typical par- 
ent species ejected from grain mantles. In the log-log figures 
presented here the relative HNCO abundance is seen to linearly 
increase with time (the actual abundance is increasing exponen- 
tially with time), reaching a peak value of ~ 10"^, consistent 
with the peak observed values, after about 2x10^ years (Fig.|6l). 
This growth suggests that the gas phase HNCO is in fact a daugh- 
ter product from the species ejected from the grain mantles. 

Examining the grain mantle abundances shows that although 
it is efficiently formed, HNCO is also rapidly processed to more 
complex species. It is the destruction of these species after their 
ejection which produces the gas phase HNCO. In our models the 
HNCO is allowed to form several daughter products (HNCHO, 
HNCOCHO, HNCONH and HNCOOH) on grain surfaces dur- 
ing the collapse phase. Once the hot core stage is reached it 
is these species that are ejected and subsequently destroyed in 
the gas-phase thus retur ning HNCO. Such a pos sibility has been 
previously suggested bv lZinchenko et all (|200C|) , although these 
authors prefered the explanation that the HNCO they observed 
was formed in the gas phase in post-shock gas. If the HNCO 
in our models was not destroyed on the grains it would be re- 
turned to the gas-phase as soon as the temperature is increased 
to 200 K with a relative abundance of ~ 2 x 10"^'. This abun- 
dance is several orders of magnitude overabundant compared to 
observations. 

Based on its excitation temperature iBisschop et al.l (l2007h 
argue that HNCO is a 'first generation' species ejected from 
grains. Our models indicate that this is not the case and sug- 
gest that this association of a molecular species with physical 
conditions on the basis of their excitation temperature alone can 
be misleading. 

4.1. HNCO/CS abundance ratio 

The HNCO/CS abundance ratio as a function of time is shown 
for the PDR models along with all the HC models in Figure |9] 
The PDR models never reach values greater of the HNCO/CS 
abundance ratio > 10 orders of magnitude smaller than ob- 
served. For model HCl the ratio increases from ab out 2 x 10'** 
at 10^ years, peaking at ~ 3 around 3 x 10^ years. iMartm et al] 
(120081) observed ratios for HNCO/'^CS in the range ~ 1 to 
~ 100. Adopting a valu e of 77 for the abundance ratio to C/'^'C 
(IWilson&Roodlll994l) . this would imply HNCO/CS ratios of 
~ 0.013 to ~ 1.3. The models span this ratio for times between 
approximately lO'* years to 2 x 10^' years showing that the time 
dependent evolution of the chemistry of hot cores alone can pro- 
duce the entire range of observed HNCO/CS abundance ratios. 
Therefore, independent of its effect on the ratio, it is clear that 
HNCO/CS is not a unique tracer of UV radiation. 

As the figure also shows, the underlaying sulphur abundance 
also has a direct effect on the absolute value of the HNCO to CS 
abundance ratio. When S is reduced by a factor of 10, for a given 
time less than about few times 10^ years, the HNCO/CS ratio 
increases by a factor of 10. However for these times the evolution 
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Fig. 9. HNCO/CS abundance ratio as a function of time for PDR 
models (top), hot core models 1-4 (middle) and hot core models 
5-10 (bottom). The shaded region defin es the upper and low er 
limit for the HNCO/CS ratio as found by lMartm etaP (l2008h . 



for different sulphur abundances is essentially independent of 
abundance. Overall the ratio monotonically increases with time 
up to a peak at a few times 10^ years. This generally mono tonic 
behaviour suggests that for a given sulphur abundance this ratio 
could act as a 'chemical clock' to constrain the ages of sources. 
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5. Summary 

Our model shows that gas phase reactions appear incapable of 
producing HNCO in sufficient quantities (even at the ~ 200K 
temperature of hot cores) to match the observed abundance of 
this species. However the rate of a key gas phase reaction is 
uncertain and new laboratory determination of this rate coeffi- 
cient is highly desirable. HNCO can be formed on icy surfaces 
by a variety of reactions. Despite the existence of these plausi- 
ble mechanisms which can produce sufficient HNCO to match 
observations, laboratory data is again missing and should be ob- 
tained. Although it is produced on the grain surfaces in poten- 
tially large amounts, our models show that the gas phase HNCO 
is in fact a daughter product from the breakdown of more com- 
plex species to which the HNCO on the grains have been pro- 
cessed. This is contrary to recent speculation about the origin of 
this species (' Bisschop et al.l 1200 7). The model also shows that 
in PDRs HNCO can not be directly produced in sufficient quan- 
tities to match observations. In hot cores in particular it seems 
that surface chemistry during an earlier cold collapse phase fol- 
lowed by gas-phase processing (in the post-collapse regime) is 
important in the production of HNCO. 

The ratio of HNCO to CS evolves strongly as a function of 
time while its absolute value is a sensitive function of the sulphur 
abundance. Whilst observations have suggested this ratio can be 
used to trace different physical conditions, our findings show that 
metallicity, choice of chemical network and hot core age can 
all also contribute to this value. Indeed, considering that HNCO 
contains three heavy elements one may expect the abundance of 
HNCO to be highly sensitive to the metallicity of a region. 
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